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Many technological applications of multiferroics are based on their ability to reconstruct the 
domain structure (DS) under the action of small external fields. In the present paper we analyze 
the different scenarios of the DS behavior in a multiferroic that shows simultaneously ferro- and 
antiferromagnetic ordering on the different systems of magnetic ions. We consider the way to control 
a composition of the DS and macroscopic properties of the sample by an appropriate field treatment. 
We found out that sensitivity of the DS to the external magnetic field and the magnetic susceptibility 
in a low-field region are determined mainly by the destressing effects (that have magnetoelastic 
origin). In a particular case of S^CuaCUCk crystal we anticipate the peculiarities of the elastic and 
magnetoelastic properties at T ~ 100 K. 
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I. INTRODUCTION 

During the last ten years a special attention is paid to the materials in which magnetism coexists with the other 
types of ordering, i.e., ferroelectric 1 , elastic 2 , martensitic . Solids that show strong coupling between the different 
types of ordering are often called multiferroics 1 . Growing interest to multiferroics is based on the possibility to i) 
control such macroscopic properties of a sample as conductivity, magnetization, elongation, with the suitable fields of 
different nature; ii) manipulate the state of the magnetically (electrically, etc.) inert materials (see, e.g., Refs. 5 and 
6). 

One of the most technologically appealing property of multiferroics, namely, sensitivity of their macroscopic prop- 
erties to the influence of small external fields is due to formation and reconstruction of the domain structure (DS) . 
This adaptivity, ability to change macroscopic parameters (such as a shape, magnetization, electric polarization) in 
response to external forces is related with the finite size and boundary of the sample. While the physical mechanism 
of the DS formation is related with the sample boundary, reconstruction and restructurization of the domains under 
external fields depends upon the properties of the domain walls. If a potential barrier for the domains wall formation 
is high, switching between the different macroscopic states is sharp and field dependence of macroscopic parameters 
reveals a hysteresis. In the opposite case of low potential barrier, reconstruction of the DS takes place through the 
nucleation and growth of new domains and shows the features of liquid-like behavior: nonhysteretic transitions be- 
tween the different macroscopic states, shape deformation, etc. The most interesting case on which we concentrate 
our attention in the present paper lies in-between: in multiferroics with the domains of different nature some types of 
the domains can easily nucleate and show soft-like behavior while the others could have high nucleation barrier and 
reveal themselves as solids. 

The origin of the DS in the "single" ferroics, like ferromagnets (FM) and ferroelectrics, is well established* and is 
attributed to the presence of long-range interactions between the magnetic or electric dipoles localized on different 
sites. The nonlocal character of the dipole forces ensures strong dependence of the equilibrium DS on the shape of 
the sample. In many important cases the DS of a single ferroic consists of the domains with opposite (sometimes 
noncollinear) directions of polarization and can be described thermodynamically on the basis of demagnetization 
energy. 

Antiferromagnets (AFM) give an example of the more complicated materials with usually pronounced coupling 
between the magnetic order parameter and lattice strain". The behavior of equilibrium domain structure in AFM is 
very similar to the behavior of the DS in other ferroics with some "technical" distinctions: i) long-range dipole-dipole 
interactions responsible for the formation of equilibrium DS have a magnetoelastic origin and are described by the 
destressing 10 (in contrast to depolarization or demagnetization) energy; ii) DS consists of the domains with different 
(nonparallel) orientations of AFM vectors. 

Description of the DS in multiferroics seems to be much more complicated problem mostly due to the fact that the 
domains of different nature appear at different scales and thus form a hierarchical structure. Good example of such 
a complexity is given by FM martensites (see, e.g., Ref. 3) that are usually characterized by two independent order 
parameters 11 , magnetization and spontaneous strain 12 . Cross-correlations between the magnetic and structural order 
parameters open a way to control a martensitic DS and, as a consequence, macroscopic deformation of a sample, with 
the external magnetic field (so-called giant magnetostriction 13 ). 

Another example of multiferroic behavior is given by some of the high-temperature superconducting systems (like 
Sr 2 Cu304Cl2 or Ba2Cu304Cl2) that show simultaneously FM and AFM ordering on the different systems of copper 
ions. In contrast to FM martensites, the DS in these crystals is not hierarchical. Each type of domains is characterized 
by two independent (FM and AFM) order parameters. Though macroscopic state of both ferro- and antiferromagnets 
could be controlled by the same, magnetic, field, the responses of the FM and AFM domain structures are different, 
as illustrated in Fig. 1. The domain structure of FMs reconfigures in the magnetic field which is parallel to an easy 
axis and does not change if the magnetic field is perpendicular to this axis. Macroscopic magnetization of the sample 
(and hence, macroscopic susceptibility) is inversely proportional to the appropriate component of demagnetization 
tensor. In contrast, in antiferromagnetic crystals the DS reconfigures for both mutually perpendicular orientations 
of the magnetic field. Macroscopic magnetization depends upon the components of destressing tensor that have a 
magnetoelastic origin. So, a material that bears simultaneously the features of FM and AFM can show some new 
type of behavior in the external magnetic field governed by competition between the demagnetizing and destressing 
effects. 

In the present paper we study an equilibrium DS of multiferroic Sr2Cu304Cl2 with the FM and AFM order 
parameters. In the framework of phenomenological approach we analyze the possible magnetization curves that could 
be obtained for the samples of different shape and different field treatment. On the basis of the developed model 
we make an attempt to interpret the unusual behavior of macroscopic magnetization observed in the experiments of 
Parks et al 1 1 and predict a peculiarity of the elastic properties of Sr2Cu304Cl2 at the temperature T w 100 K. 
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Figure 1. (Color online) Behavior of the FM (a) and AFM (b) domain structure in the external magnetic field H. In the 
absence of field both types of domains (shown by arrows) are equally represented, (a) FM domains have opposite direction of 
magnetization vector. The magnetic field applied parallel to an easy axis (upper panel) removes degeneracy of the domains. 
As a result, fraction of the favorable domain increases. If H is perpendicular to the easy axis, domains of both types are 
equivalent, the domain fraction does not change, magnetic field induces a tilt of the magnetizations (lower panel), (b) AFM 
domains have different (perpendicular) orientations of AFM vectors. Degeneracy of the domains is removed for any of two 
mutually perpendicular orientations of the magnetic field. 



II. MODEL 

The crystal structure of high-temperature superconducting cuprates Sr2Cu304Cl2 and Ba2Cu304Cl2 consists of 
CU3O4 planes separated by spacer layers of SrCl or BaCl 14-11 '. Two types of magnetic ions, Cul and Cull (see Fig. 2) 
form two interpenetrating square lattices within CU3O4 planes. Within the temperature interval Tjj = 40 K< T < 
Ti = 380 K the ions of the first type (Cul) are AFM ordered while the ions of the second type (Cull) bear small but 
nonzero FM moment . According to the experiments 18, mutual orientation of Cul and Cull moments depends upon 
the direction of the external magnetic field and can be either perpendicular or parallel. Thus, the magnetic structure 
consists of two weakly coupled subsystems, namely, an AFM, localized on Cul ions, and a FM one, localized on Cull 
ions. The FM subsystem is unambiguously described by the magnetization vector Mp and the AFM subsystem is 
described by two vectors: AFM vector L = (Si — S2 + S3 — S 4 )/4 and ferromagnetic vector M = J2j Sj/4 (numeration 
of Cul sites is shown in Fig. 2). 

In the absence of the external field the FM moments at Cull sites are oriented along (110) crystal directions 
perpendicular to the staggered magnetizations of AFM subsystem, as shown in Fig. 2. Due to tetragonal symmetry 
of the crystal (space group IA/mmm) an equilibrium magnetic structure can be realized in four types of equivalent 
domains as shown in Figs. 2 and 3. Domains of type A and B could be thought of as AFM domains because they 
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Figure 2. (Color online) Magnetic structure of CU3O4 layer in two different configurations (domains). Magnetic field is parallel 
to (110). Two types of magnetic ions are represented with the filled and hollow circles. FM ordered moments of Cull could be 
(a) parallel (domain A) or (b) perpendicular (domain B) to the applied magnetic field. Small canting of the Cul spins induced 
by the external magnetic field is not shown. 



correspond to different orientations of L vector and thus are sensitive to orientation of the magnetic field H with 
respect to the crystal axes (see Fig. 1). Types Al and A2 (and, correspondingly, Bl and B2) are FM domains, they 
have an opposite direction of Mp vector and could be removed from the sample by H||Mp. 




I [HO] 



Domain Bl, p=3n/4, 9=nl4 



Domain B2, <p=-n/4, 9=5k/4 



Figure 3. (Color online) Four types of magnetic domains. Axes x and y are parallel to (100) crystal directions. The external 
magnetic field H||[110] (if any). Types A and B have different orientations of AFM vector, types 1 and 2 correspond to opposite 
directions of FM vector M_p. Ellipse (dash line) images the supposed shape of the sample and its orientation (axes X, Y) with 
respect to crystal axes. 

Phcnomenological description of the DS is based on the analysis of the free energy potential $ of the sample. We 
take into account three constituents of $: magnetic $ ma g, stray (demagnetizing), $ stray and destressing, $dost, energy: 

* = $mag + Astray + 3>dcst- (1) 

Magnetic energy of Sr2Cu304Cl2 crystal in mean field approximation is well established 10 ' 18,19 and can be written as 
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Table I. Parameters used in the free energy [Eq. (2)]. The second column gives the raw data (in meV) as taken from Refs.14, 
18, and 19, in the last column the same values are given in Oe. 



Parameter 


Meaning 


Value in meV 


Value in Oe 


Jo 


Cul-Cul superexchange (in-plane) 


130 


1.02-10 7 


J av 


isotropic pseudodipolar interaction 


-12 


-9.4-10 5 


Jpd 


anisotropic pseudodipolar interaction 


-0.027 


-2.1-10 3 


K± 


out-of-plane anisotropy 16 


0.068 


5.3 • 10 3 




in-plane anisotropy 


10-10 -6 


7.8-10" 2 



follows: 



$ mag = f dV f-ij [J (M 2 - L 2 ) + J av M ■ M F + J pd M F ,T 2 L + K^L 2 ] 
Jv I M o 



Ml 



K\\L 2 x L 2 y - H • M F - 2H • m| 



(2) 



Here V is the sample volume, Mo is Cul sublattice magnetization, orthogonal axes x and y are parallel to the 
crystal directions [100] and [010], respectively (see Fig. 3). a z is the Pauli matrix. The meaning and values of 
phenomenological constants are given in Table I. In the last column of this Table all the constants are converted to 
Oe by division by sublattice magnetization Mo = 27.4 Gs (that corresponds to spin s = 1/2 per Cul site). 

Contributions $ stray and $dcst in Exp. (1) arise from the long-range dipole-dipole interactions of the magnetic and 
magnetoelastic nature, correspondingly, and depend upon the sample shape. We consider a thin (thickness c) pillar 
with an elliptic crossection whose principal axes X and Y are parallel to (110) directions within the CU3O4 layers 
(see Fig. 3). In this case 

Astray = j [iV a dm (M FX + 2M X ) 2 + N^{M FY + 2M Y ) 2 ] , (3) 

where the brackets (. . .) mean averaging over the sample volume. The components of demagnetization tensor iV^™ 
are calculated in a standard way : 

Ndm = 4ttc r' 2 sin 2 ^ dm = AncVT^ cos 2 g# 

aVT^Jo y/l-k 2 sinV " « Jo ^l-k 2 '/ " ' ' 



■ sm 1 



Here a > b(^> c) are the ellipse's semiaxes (parallel to X and Y axes) and the parameter k 2 = 1 — b 2 /a 2 depends 
upon an aspect ratio b/a of the sample. 

The destressing energy can be written in an analogous form 

$dcs t _ v {^[(if, - L\) 2 + 4{L x L- x2 



Y 

+ Nii«(L x L 2 Y ) <-[(L| - L 2 Y ) 2 A{L X L Y ) 2 ]} . (5) 

An explicit form of the destressing constants 7V dcs depends upon the elastic and magnetoelastic properties of the 
crystal which we assume to be isotropic (that means, in particular, the following relation between the elastic modula: 
en - C12 = 2C44). Then, 



dcs _ A 2 (3-4zv) /V des = £ [A 2 (2-3^) + A,A]J 2 (fc) dcs c A 2 J 4 (fc) 
is 16cm(1-i/)' 2 b 8c44(l-i/) ' 4an b 3c 44 {l - v)' 



where A and A,; are magnetoelastic constants,^ = cyijicw + c\%) is the Poisson ratio and we have introduced the 
dimensionless shape- factors Ja,4(fc) as follows 10 



Mk) 



r/2 (sin 2 <? !) + cos2(/.//fc 2 )d0 
\J\ — k 2 sin 2 <f> 



. T/2 (1-8 cos 20 - k 2 sin 2 <j) + 8 cos 20/fc 2 )rf0 
VI - k 2 sm 2 <j) 
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In Eqs.(3) and (5) we have omitted ^(||z| [001])-components of the demagnetizing and destressing tensors as inessential 
for further consideration. 

Expressions (2), (3), and (5) could be substantially simplified if one takes into account that: i) far below the 
Neel temperature the values of sublattice magnetizations Mq and M-p are saturated and constant; as a result ii) 
L _L M and L 2 + M 2 = M§ (normalization conditions); Hi) if the magnetic field is much smaller than the spin- 
flip field, H <C Jo /Mo and coupling between the FM and AFM subsystems is much smaller than AFM exchange, 
J av Mp <C JoM), the magnetization induced in AFM subsystem is small, M <C Mo, and vector M can be excluded 
from Eq. (2) as follows 21 : 



L x 



(8) 



iv) if out-of-plane anisotropy is strong enough, K± 3> K\\ (see Table I), all the magnetic vectors lie within xy (and, 
equivalently, XY) plane and could be described with the only angle variable, as shown in Fig. 3: 

L x — MocosQ, L y = Mo sin 9; Mp x = mpMoCOSip, Mp y = mpMosmip. (9) 

Here tof(=10 -3 for Sr2Cu304Cl2 ) is a dimensionless constant that represents the ratio between the spin moments 
localized on Cull and Cul sites. 

With account of the relations (8) and (9) the specific potential <f> = &/V (see Exp. (1) takes the following form 

J 2 

(j) = 4J pd TO F (cos(6» + if)) +i^||(cos46») - -^mp(cos2(6> - ip)) 

oJr 







m-pH 



1 - £L\ (cosfa - $)) + ||<cos(20 - $ - <p)) 



+ 32^0 <COS 2{6 _V0> 



- <-(cos 2(9 - + ±M m F [A dm (cos(^ - </>)) 2 + A^sin^ - V>)} 2 ] (10) 

+ N dcs (cos2(9 - V)) 2 + A7V dcs (sin 2(9 - ?/>)) 2 . 

where V is an angle between the magnetic field and x-axis, iV dcs = iV dcs + A^ 4 d a c ^, AiV dcs = 4 (iV dcs - iV 4 d a c ^), and 
we assume that the field is parallel to one of the principal axes of the sample (this corresponds to the experimental 
situation that will be discussed below). Here and for the rest of the paper we use the values in Oe (see the last column 
of Table I) instead of energy units (say, cf) — > <j>/M , etc.). 

Let us consider the case when the magnetic field is parallel to one of the easy axes, H||[110], so, ip — 7r/4. In an 
infinite sample (all the components of tensors iV dm , N des are equal to zero) minimization of 4> with respect to magnetic 
variables 9 and ip gives rise to the four solutions labeled as Al,2 and Bl,2 (see Fig. 3). Equilibrium values at H = 
are 



state Al 



state Bl 
state A2 
state B2 



= -7r/4, <pai = tt/4; 



0bi = tt/4, <pai = 3?r/4; (11) 
6»A2 = 37r/4, ip A2 = 57r/4; 

9b2 = 57T/4, (fB2 = -7t/4. 

It should be stressed that in contrast to pure AFMs the configurations with (Mir, L) and (Mf, — L) are inequivalent, 
due to anisotropic pseudodipolar interactions (described by the constant J p d). 

Fig. 4 illustrates the field-induced variation of equilibrium magnetic configurations (represented by A-projections 
of Mf and L vectors) obtained from the numerical minimization of the Exp. (10) using the data from Table I. It is 
clearly seen that within the interval \H\ < i? s _fi = 525 Oe there exist all four states Al,2 and Bl,2. The magnetic 
field removes degeneracy between the states Al, A2 and B 22 , as can be seen from Fig. 4a. In particular, when H > 0, 
the specific energies (f>j = 4>(9j,tpj) of equilibrium states are related as follows: 0ai < 4>b < 4>A2- So, in some cases 
(discussed below) variation of the external field may induce formation of the AFM (B) instead of the FM (A2) domain. 
Orientations of Mp and L vectors in the A states are not influenced by the field, while in the B states both vectors 
are slightly tilted (see Fig. 4b,c). Rotation of AFM vector from the field direction in the state B (where H||L) is a 
peculiar feature of the FM+AFM multiferroic caused by pseudodipolar interactions between Cul and Cull ions. In 
the pure antifcrromagncts an AFM vector L keeps parallel (with respect to H) orientation up to the field of spin-flop 
transition. 

The first critical field i/ s -fi ex JoK\\ corresponds to a step-like (spin-flop) transition Bl,2— >Al. In the interval 
i/s_fi < \H\ < H s -{2 = 1465 Oe the potential $ has only two minima that correspond to the states Al and A2. The 
second critical field H s _f 2 corresponds to 180° switching of Mp vector. Its value depends on the effective anisotropy 
that originates from the pseudodipolar coupling (corresponding constants J aV j^pd) and in-plane anisotropy K\\ and 
can be calculated only numerically. Above H > i? s -f2 the sample is in a single domain state (Al). 
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Figure 4. (Color online) Stability ranges of homogeneous configurations shown in Fig. 3 in the external magnetic field H||[110]. 
(a) Specific energy (in Oe) of equilibrium homogeneous state vs H. (b), (c) Normalized projections of FM and AFM moments 
on the field direction. Field induces rather noticeable rotation of Mp vector toward the field direction (b) and slight tilt of L 
vector (c). Loss of stability takes place at the critical values H = Hs-np, as shown with arrows. 



III. EQUILIBRIUM DOMAIN STRUCTURE AND MAGNETIZATION CURVES 



On the large scales (much greater than the characteristic scale of the magnetic inhomogeneity, i.e., domain wall 
thickness) the magnetic structure of the sample is represented by a set of magnetic variables {6j, ipj} that describes 
orientation of FM and AFM vectors inside the domains (j =A1, A2, Bl, B2) and a set of variables {£j} that represents 
the amount of matter (say, volume fraction) in the state of j-type (obviously, J^C? = !)■ Equilibrium DS in presence 
of the external field is then found from the condition of minimum of $ with respect to {£?_,-, <fj, 
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In such an approach one can neglect a contribution of the domain walls into free energy potential $. However, 
we implicitly account for the inhomogeneities in space distribution of the FM and AFM vectors when we chose 
independent variables for the potential $. Namely, reconstruction of the DS may proceed in two ways: i) by the 
field-induced motion of the domain walls; ii) by nucleation and growth of the energetically favorable domains. The 
first way is almost activation-less while in the second case the system should overcome the potential barrier related 
with the formation of the domain walls. In the case under consideration the domain walls between AFM (A/B) and 
FM (A1/A2 or B1/B2) domains have different energies, and so, appear at different conditions. In what follows we 
consider some typical situations and show the way to control the DS with appropriate treatment of the sample. 



A. Four types of domains 



In the case when all four types of domains may freely grow or diminish in size (say, in a virgin sample that initially 
contains domains of all types), the external magnetic field is screened by an appropriate domain configuration (see 
Fig. 5) and the effective field inside the sample is zero. Equilibrium values of the magnetic variables in this case are 
given by Eq. (11) and the domain fractions depend on magnetic field as follows: 



6u,A2 — - 



1 - e (0) ± 



2H 



( H 



\Hdcs 



£B1 - S.B2 — J 



(0) _ 



H 

Hdcs 



where we have introduced the demagnetizing, ffj m , and the destressing, Hd cs , fields: 



H dm = m F 7V Q dm M , H dcs = 8y/JoN des . 



(12) 



(13) 



As seen from Eq. (13), the value of destressing field is enhanced due to exchange interactions (constant Jo). On 
the contrary, the demagnetizing field is weakened due to small FM moment (mp <C 1). So, in the crystal under 
consideration the demagnetizing effects are much smaller than the destressing ones, -f/dm <C -ffdcs (see Table III A). 

The value £W introduced in Eq. (12) represents the disbalance between type A and type B domains in the absence 
of field. This value depends upon the shape of the sample (or, equivalently, from the aspect ratio, see Eq. (6)): 



Aides 
(0) _ JV 2an 



- -Ja(fc). 



(14) 



Such a shape- induced nonequivalence of domains has a magnetoelastic origin (see Ref.10 for details) and originates 
from the AFM properties of the system. The disbalance between type A and type B domains was noticed in the 
experiments Ref.14 for the different sample shapes. The value = 0.22 calculated from Eq.(14) for the typical 
sample size (see Table III A) fits well the experimental magnetization curves, as we will see below. 

The described configuration of the DS (see Eq. (12)) is schematically shown in Fig. 5b. The fraction of the 
unfavourable domains A2, Bl and B2 diminishes and at the critical value 



H = H C 



H, 



dm 



i-^/i-(i-£ (0) )( 



Hdr, 



r(l-£ (0) )ffdr 



(15) 



the unfavourable FM domains A2 disappears (£a2 = 0). 

At H > H cl i the internal effective magnetic field is nonzero and magnetizations in the domains of B type rotate. 
However, if iJ cr i -C i? s _f (as, indeed the case in the crystal under consideration), small tilt of Mp and L vectors can 
be neglected and field dependence of the domain fractions (shown in Fig. 5c) is approximated as 



Al 



(0) 



16-Htof Jo 

H dcs 



H 

Hdc: 



Bl,2 



(0) 



16-ffmp Jo 

H dcs 



H 

Hdc: 



(16) 



The second critical field at which the unfavourable domains of B type disappear (£bi.2 

2A^ dl 



H CI 2 = 8tof Jo 



(1 + (m) - 1 



771F 



0) is given by the expression 

(i + e (0) ). (17) 



Above the second critical field, H > H cr2 , the sample is a single domain (Al) in average, with the possible remnants 
of the states A2, B and corresponding domain walls that can serve as the nucleation centers during the field cycling. 
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Figure 5. (Color online) Behavior of the combined FM and AFM domain structure in the external magnetic field H||[110] 
(parallel to the long side of the sample), (a) In the field absence all types of domains (Al, A2, Bl, B2) are equally represented, the 
disbalance £' ' between A and B types depends upon the aspect ratio of the sample, (b) Below the first critical field, H < Hal, 
the domains are rearranged in such a way that the effective magnetic field vanishes, (c) In the interval H CI \ < H < H CT 2 the 
unfavourable domains of A2 type disappear, Al-domains compete with the domains of B-type. (d) Above the second critical 
field, H > H cr 2, the sample is a single domain. 



Table II. Parameters used 


in numerical simulations. The source of data 


(experimental or calculated) 


is specified in the last 


column. 








Parameter 


Meaning 


Value 


Rem 


a x b x c 


Sample size 


7 x 2 x 0.5 mm 3 


Ref.14 


niF 


Mf/Mo 


7-10" 4 


Ref. 


M F 


Saturation magnetization 


7 • 10~ 3 emu/g 


Ref. l I 


£(0) 


Shape-induced bias 


0.22 


Eq.(14) 




Demagnetization field 


0.3 Oe 


Eq.(13) 


N dcs 


Dcstressing const., T = 120K 


7 mOe 


Fitting 




T = 100 K 


1.5 mOe 


par am. 


Hdcs 


Destressing field, T = 120 K 


2.1 kOe 


Eg. (13) 




T — 100 K 


1.1 kOe 





Full monodomainization of the sample takes place above the critical field -ff s -f2 3> H CI 2 at which all the states except 
Al became unstable. 

Field cycling of the sample that initially had all the types of domains is reversible if the maximal field value H max 
is not very large, H CT 2 < H m&x <§C -ff s _f2- Macroscopic magnetization is parallel to the direction of the external field 
due to the full compensation of the perpendicular component by Bl and B2 domains. 

Field dependence of macroscopic magnetization Af par oc (£ai — £a2) at T — 120 K calculated from Eqs. (12) and 
(16) (see Tables I and III A) is represented in Fig. 6. One can distinguish three intervals that correspond to different 
domain composition: i) steep growth of M par from to oc 0.5(1 — ^")Jlfp (at H = H CI i) due to the motion of A1/A2 
domain walls initiated by demagnetization; ii) smooth growth of M par from oc 0.5(1-£ (0) )Mf to w Mp (at H = H cr2 ) 
due to the motion of Al/B domain walls initiated by the destressing; Hi) very smooth growth of M par due to rotation 
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of sublattice magnetizations (not seen in Fig. 6). Such a behavior contrasts with a "standard" magnetization curve 
of FM and also with the case when only two types of domains could compete under the action of external field. The 
last case will be considered in details in the next section. 



A1+B 




H(Oe) 



-1.0 -J 



Figure 6. (Color online) Magnetization curve (projection on field direction) in the external magnetic field H||[110] for the case 
depicted in Fig. 5. Inset shows the details of magnetization behavior below H CT i. The jogs (shown with arrows) arise at the 
critical fields H — H C ri,2 when one type of domains disappear. Magnetization is normalized to saturation value. 



B. Competition of two domains 

Let us consider a sample that was preliminary monodomainized to the state Al by excursion into the region of high 
field, H > i? s -f2- The DS in this case depends upon the relation between nucleation energies of different states. As it 
was shown above, at H > an AFM domain B is more favorable than a FM domain A2 (see Fig. 4a). If, in addition, 
there is a slight misalignment between the magnetic field H and a crystal axis [110] that removes degeneracy between 
the Bl and B2 states, the DS of the sample is represented by the domains of only two types, Al and Bl. 

Equilibrium values of the magnetic variables in this case were calculated by the numerical minimization of the 
potential (10) with limitations ^A2 = £b2 = 0. The values of the destressing coefficient N dcs at different temperatures 
(see Table III A) were defined from the fitting of experimental data 14 . 

Field dependence of macroscopic the magnetization at T = 120 K is shown in Fig. 7 with solid lines. Points 
represent experimental data 1 '. Due to the fact that the domains Al and Bl could not screen the external field, the 
macroscopic magnetization has two components: one that is parallel to H (upper panel in Fig. 7) and one that is 
perpendicular to H (lower panel). The parallel and perpendicular components represent the fractions of Al and Bl 
domains, respectively. 

When the field decreases from high positive values, an AFM domains of B type appear and magnetizations M par (i?) 
and M peTp (H ) vary smoothly between zero and saturation value. The slope of magnetization curves depends upon the 
destressing coefficient and is thus much smaller than the initial steep slope in the 4-domain case (see Fig. 6). At small 
negative field the sample is almost a single domain (type B). However, this state is a metastable one from the energy 
point of view, as seen from Figs. 8 and 4a. Really, below H = H3 w 2.8 Oe (marked with arrow in Fig. 8) the energy 
of the state A2 (with M_p tl H) is lower than that of a single domain state Bl and a multidomain state Al+Bl. On 
the other hand, due to preliminary high-field treatment, the sample contains no nucleation centers of A2 state. So, 
the states Bl and A2 are separated with the potential barrier that could be overcome only at H — H2 ~ —50 Oc 
(according to Ref. I, this value varies from sample to sample and depends on temperature). After the subsequent 
excursion into high negative fields (well below H2) the sample transforms into a single domain A2 and one can observe 
competition between the A2 and B2 domains during the further field increase. 
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Figure 7. (Color online) Macroscopic magnetization us magnetic field. Points - experimental data for Sr2Cu304Cl2 taken 
at T —120 K after monodomainization of the sample at high fields H tx 5 T. Solid lines - theoretical approximation (see text 
for details). Upper and lower panels show, correspondingly, the parallel and perpendicular components of magnetization with 
respect to magnetic field. Insets show geometry of the experiment: orientation of the field with respect to the sample and 
orientation of the measured magnetization with respect to H. The dominant type of domains for each field interval is depicted 
schematically by the single- and double-headed arrows. 



C. Domain structure and field treatment 



In the previous subsections we have considered two limiting cases of field treatment that result in two types of 
magnetization curves. In the virgin sample (no field treatment) the magnetization can be smoothly and reversibly 
changed between two opposite directions. Field cycling between high fields (high enough to remove all the domain 
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Figure 8. (Color online) Magnetic energy vs magnetic field, T = 120 K. Lines 1 and 3 correspond to a single domain state 
(domains Al and A2, correspondingly), line 2 represents equilibrium two-domain state (domains Al and Bl). Domain Bl 
appears at H — Hi. Domain A2 appears at field H = H2 (determined empirically) when the energy difference between 
two-domain and single domain states is enough for nucleation of this energetically favorable domain. At H — H3 the energy 
of two-domain state is equal to the energy of a single-domain state A2, but the potential barrier between two states prevents 
nucleation of the domain A2. 



walls and the remains of unfavorable domains) results in a hysteretic behavior when magnetization varies smoothly 
between zero and saturation value and then suddenly changes due to transition from metastable to stable state. 

In this subsection we consider some intermediate case when a single domain sample is cycled in low fields. Corre- 
sponding magnetization curve is shown in Fig. 9a (solid lines - numerical simulations, points - experimental data 14 , 
T = 100 K). Field cycling starts at high positive fields, where the sample is a single domain. When the field is 
decreased down to H = Hi (see Fig. 8a) the domains of Bl type appear and the DS consists of Al and Bl domains. 
Magnetization curve (upper curve in Fig. 9a) in this case is of two-domain type discussed in Subsection IIIB (we 
still assume slight misalignment that excludes one type of B domains). However, further behavior of the DS and 
hence, magnetization, depends upon the size of the loop. If the loop is small (\H\ < |i?a|) where H2 is a coercive 
field at which domain Bl transforms into A2 as explained above), magnetization varies smoothly between zero value 
at negative fields and saturation value at positive fields. If the loop is very large (\H\ > \H cl2 \), the DS structure 
consists of two domains: Al and Bl for large positive and small negative fields and A2 and B2 for large negative and 
small positive fields, as shown in Figs. 9b and 7. In the intermediate case (|-Ha| < \H\ <C | i? C r2 1 ) the DS includes three 
types of domains, Al, A2 and B2 (lower curve in Fig. 9a), and magnetization curve is asymmetric. It is worth to note 
that theoretical magnetization curves calculated with only one fitting parameter (destressing coefficient N des ) fit well 
experimental data, as seen from Figs. 9 and 7. 



IV. DISCUSSION 

We have considered the different types of the DS behavior in a multiferroic with AFM and FM order parameters. 
Depending on the field treatment the DS may include from one to four types of domains and can be unambigu- 
ously determined from the magnetization curves in the field H||[110]. Namely, if DS includes all types of domains, 
macroscopic magnetization is parallel to H, magnetization curve is reversible, varies between ± saturation value, and 
includes steep section at small fields. If the DS includes three types of domains Al, A2 and Bl, the macroscopic 
magnetization has two components, parallel, M paT , and perpendicular, Af pcrp , to H. During field cycling M par varies 
between positive and negative saturation values, while M perp varies between zero and saturation value. At last, if the 
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DS includes only two domains, A and B, both M par and M pcrp vary between zero and saturation value. 

We argue that formation of AFM (B) domain results from the destressing effect which, in turn, originates from 
magnetoelastic interactions. An absolute value of magnetoelastic constant is rather small (compared to such AFMs 
as NiO, KC0F3, etc) and corresponds to spontaneous strain u oc 10 -6 (for estimation we took C44 = 20 GPa at 
T = 120 K). Such a small value of u explains the low potential barrier for formation of AFM domains. 

Analysis of magnetization curves shows that low- field susceptibility \ of the sample that consists of AFM domains 
is inversely proportional to the destressing coefficient _/V dos (in contrast to FM, where x depends upon demagneti- 
zation constant). According to the experiments 14 , the inverse susceptibility x _1 of Sr2Cu304Cl2 shows nontrivial 
temperature dependence (see Fig. 10b) and attains the minimum at T = Tq — 97 K. The domain fraction £ai at fixed 
H extracted from the neutron scattering experiments of the reminded group ' shows the same temperature depen- 
dence as x, as can be seen from Fig. 10b. Using correlation between iV dcs and \ we predict the following temperature 
dependence of the destressing coefficient depicted in Fig. 10a: 

jydesfjn _ I 7.3 • 10~ 5 • (T - T ), T > T , 
\ 6.13-10- 5 -(T -T), T<T . 

If we take into account that iV des cx A 2 /c44 (see Eq.(6)) we may also anticipate a peculiarity of the elastic (or 
magnetoelastic) properties of the crystal in the vicinity of T = To. 

In summary, we described the possible scenario of field-induced restructurization of the domains in the system 
that consists of the domains of different physical nature. The proposed model can be extended to multiferroics that 
show simultaneously ferroelectric and AFM ordering and also to FM martensites with ferroelastc and ferrpmagnetic 
ordering. 
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Figure 9. (Color online) Macroscopic magnetization vs magnetic field for different field treatment. Points - experimental 
data for S^CuaCMCh; ' 1 taken at T =100 K (see text for details), solid lines - theoretical approximations. Thin arrows 
show the direction of field sweeping, (a) Competition between FM and AFM domain structure. Field treatment starts from 
the high positive values (H « 1 T) where the sample is a single domain. Variation of field is swept at intermediate value 
\H2\ < \H\ <g |-?/ cr 2|, large enough to induce switching between metastable, Bl, and stable, A2, states and small enough to 
remove the traces of Al and Bl phases from the sample, (b) Hysteresis loop with the excursion into high fields. For any field 
value the DS includes only two types of domains, as in Fig. 7. Only parallel component of magnetization is shown. 
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Figure 10. (Color online) Temperature dependence (a) of the destressing coefficient iV dcs predicted from the comparison with 
the (b) temperature dependence of the reciprocal domain Al fraction (2£ai — 1) _1 (triangles) and the inverse susceptibility x _1 
(diamonds) plotted according to data 14 . Raw data for x _1 were normalized (multiplied by appropriate factor) to fall into the 
same ranges of values as (2£ai — 1) _ • Dash line shows linear approximation of the experimental data. Peculiarity at T = 97 K 
is indicated with arrow. 



